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O ■ We study quantum effects on moduli dynamics arising from the production of particles 

^ ■ which are light at special points in moduli space. The resulting forces trap the moduli 
CL|' at these points, which often exhibit enhanced symmetry. Moduli trapping occurs in time- 
^ . dependent quantum field theory, as well as in systems of moving D-branes, where it leads 
■ the branes to combine into stacks. Trapping also occurs in an expanding universe, though 
^ I the range over which the moduli can roll is limited by Hubble friction. We observe that a 
5^ . scalar field trapped on a steep potential can induce a stage of acceleration of the universe, 
which we call trapped inflation. Moduli trapping ameliorates the cosmological moduli 
problem and may affect vacuum selection. In particular, rolling moduli are most powerfully 
attracted to the points with the largest number of light particles, which are often the points 
of greatest symmetry. Given suitable assumptions about the dynamics of the very early 
universe, this effect might help to explain why among the plethora of possible vacuum 
states of string theory, we appear to live in one with a large number of light particles and 
(spontaneously broken) symmetries. In other words, some of the surprising properties of 
our world might arise not through pure chance or miraculous cancellations, but through a 
natural selection mechanism during dynamical evolution. 
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1. Introduction 

1.1. Moduli Trapping Near Enhanced Symmetry Points 

Supersymmetric string and field theories typically contain a number of light scalar 
fields, or moduli, which describe low-energy deformations of the system. If the kinetic 
energy of these fields is large compared to their potential energy then the classical dynamics 
of the moduli is described by geodesic motion on moduli space. 
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At certain special points (or subspaces) of moduli space, new degrees of freedom 
become light and can affect the dynamics of moduli in a significant way |1|,0,^J|,0. These 
extra species often contribute to an enhanced symmetry at the special point. We will refer 
to any points where new species become light as ESPs, which stands for extra species 
points, and also, when applicable, for enhanced symmetry points. 

A canonical example is a system of two parallel D-branes. When the branes coincide, 
the two individual ^7(1) gauge symmetries are enhanced to a t/(2) symmetry, as the strings 
that stretch between the branes become massless . Similar points with new light species 
arise in many contexts; examples include the Seiberg-Witten massless monopole and dyon 
points in A/" = 2 supersymmetric field theories , the conifold point and ADE singu- 
larities in Calabi-Yau compactification 0, the self-dual radius of string compactifications 
on a torus, small instantons in heterotic string theory and many other configurations 
with less symmetry. 

Classically, there is no sense in which these ESPs are dynamically preferred over other 
metastable vacuum states of the system. We will argue that this changes once quantum 
effects are included. In particular, quantum particle production of the light fields alters 
the dynamics in such a way as to drive the moduli towards the ESPs and trap them there. 

The basic mechanism of this trapping effect is quite simple. Consider a modulus (j) 
moving through moduli space near an ESP associated to a new light field x- For example, 
(f) could be the separation between a pair of parallel D-branes, and x ^ string stretching 
between the two branes - in this case the ESP (p = is the point where the branes coincide 
and X becomes massless. As (p rolls through moduli space, the mass of x changes; x S^ts 
lighter as (p moves closer to the ESP and heavier as (p moves farther away. This changing 
mass leads to quantum production of x particles; as (p moves past the ESP some of its 
kinetic energy will be dumped into x particles. As (p rolls away from the ESP, more 
and more of its energy will be drained into the x sector as the x mass increases, until 
eventually (p stops rolling. At this point the moduli space approximation for (p has broken 
down, and all of the original kinetic energy contained in the coherent motion of (p has 
been transferred into x particles, and ultimately into all of the fields interacting with x 
(including decoherent quanta 5(p). As we will see in detail, the x excitations generate a 
classical potential for (p which drives the modulus back toward the ESP and traps it there.S 



^ There are also corrections to the effective action for (p from loops of x particles, including both 
kinetic corrections and a Coleman- Weinberg effective potential. Both effects will be subdominant 
in the weakly-coupled, supersymmetric, kinetic-energy dominated regimes we will consider. 
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In the example of the pair of moving D-branes, the consequences of this are simple: 
two parallel branes that are sent towards each other will collide and remain bound to- 
gether. The original kinetic energy of the moving branes will be transferred into open 
string excitations on the branes and eventually into closed string radiation in the bulk. 

In §2 we will describe the general trapping mechanism and study its range of appli- 
cability using a few simple estimates. In §3 we will write down the equations of motion 
governing trapping in more detail, and describe the numerical and analytic solutions of 
these equations in a variety of cases. 

It is important to recognize that this trapping effect is in no way special to string 
theory. Flat space quantum field theory with a moduli space for (f) and an ESP is an ideal 
setting for the trapping effect, and it is in this setting that we will perform the analysis of 
§2 and §3. In §4 we will generalize this to incorporate the effects of cosmological expansion, 
and in §5 we will discuss the possibility of significant effects from string theory. Having 
established the moduli trapping effect in a variety of contexts, we will then study its 
applications to problems in cosmology. 

The most immediate application is to the problem of vacuum selection. As we will 
see in §6, the trapping effect can provide a dynamical vacuum selection principle, reducing 
the problem to that of selecting one point within the class of ESPs. This represents 
significant progress, since the vast majority of metastable vacua are not ESPs. Trapping 
at ESPs may also help solve the cosmological moduli problem, as we will see in §7. In 
particular, trapping strengthens the proposal of by providing a dynamical mechanism 
which explains why moduli sit at points of enhanced symmetry. 

Finally, as we will explain in §8, the trapping of a scalar field with a potential can 
lead to a period of accelerated expansion, in a manner reminiscent of thermal infiation 
||12|| . This effect, which we will call trapped infiation, can occur in a steeper potential than 
normally admits such behavior. 

From a more general perspective, moduli trapping gives us insight into the celebrated 
question of why the world is so symmetric. The initial puzzle is that although highly 
symmetric theories are aesthetically appealing and theoretically tractable, they are also 
very special and hence, in an appropriate sense, rare. One expects that in a typical string 
theory vacuum, most symmetries will be strongly broken and most particles will have 
masses of order the string or Planck mass, just as in a typical vacuum one expects a large 
cosmological constant. Vacua with enhanced symmetry or light particles should comprise 
a minuscule subset of the space of all vacua. 
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Nevertheless, we observe traces of many symmetries in the properties of elementary 
particles, as spontaneously broken global and gauge invariances. Moreover, all known 
particles are hierarchically light compared to the Planck mass. Given the expectation that 
a typical vacuum contains very few approximate symmetries and very few light particles, 
it is puzzling that we see such symmetries and such particles in our world. 

For questions of this nature, moduli trapping may have considerable explanatory 
power. Specifically, the force pulling moduli toward a point of enhanced symmetry is 
proportional to the number of particles which become massless at this point, which is 
often associated with a high degree of symmetry. This means that the most attractive 
ESPs are typically the ones with the largest symmetry, and rolling moduli are most likely 
to be trapped at highly symmetric points, where many particles become massless or nearly 
massless. Moreover, the process of trapping can proceed sequentially: a modulus moving 
in a multi-dimensional moduli space can experience a sequence of trapping events, each 
of which increases the symmetry. These effects suggest that the symmetry and beauty we 
see in our world may have, at least in part, a simple dynamical explanation: beauty is 
attractive. We will discuss this possibility in §6. 

1.2. Relation to Other Works 

Similar effects have been described in the literature. There has been much work 
on multi-scalar quantum field theory in the context of inflation, especially concerning 
preheating in interacting scalar field theories. Some of our results will be based on the 
theory of particle production and preheating developed in the series of papers [0|2|,3, 
which explores many of the basic phenomena in scalar theories of the sort we will consider. 
Likewise, Chung et al. have explored the effects of particle production on the inflaton 
trajectory and on the spectrum of density perturbations. Although we will derive what we 
need here in a self-contained way, many of the technical results in this paper overlap with 
those works, as well as with standard results on particle production in time-dependent 



systems as summarized in e.g. [|T^. Although we will not study the case in which x goes 
tachyonic for some range of (f), our results may nevertheless have application to models of 
hybrid inflation p^JTSd , including models based on rapidly-oscillating interacting scalars 



0,0,0. 

In strong 't Hooft coupling regions of moduli spaces which are accessible through the 
AdS/CFT correspondence, virtual effects from the large numbers of light species dramati- 
cally slow down the motion of (f) as it approaches an ESP, with the result that the modulus 
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gets trapped there . This also provides a mechanism for slow roll inflation without very 
flat potentials. In the present work, which applies at weak 't Hooft coupling, it is quantum 
production of on-shell light particles which leads to trapping on moduli space. 

Other works in the context of string theory have explored the localization of moduli 
at ESPs. The authors of [|T9| , [20[] studied the evolution of a supersymmetric version of the 
4> — X system arising near a flop transition using an effective supergravity action. They 
showed that, given nonvanishing initial vevs for both (p and x, the flelds will settle at the 
ESP even if one formally turns off particle production effects. Our proposal, by contrast, is 
to take into account on-shell quantum effects which dynamically generate a nonzero (x^). 



In works such as |^T[] attention was focused on the boundaries of moduli space, while here 



we focus on ESPs in the interior of moduli space. In [^, production of light strings was 
studied in the context of DO-brane quantum mechanics; as we explain in §2.3, this has some 
similarities, but important differences, with our case of space-flUing branes. Scattering of 
Dp-branes was also studied in [l23|| . 



Dine has suggested that enhanced symmetry points may provide a solution to the 
moduli problem, as moduli which begin at an enhanced symmetry minimum of the quan- 
tum effective potential can consistently remain there both during and immediately after 
inflation [|ll|]. One would still like to explain why the moduli began at such a point. As 
we discuss in §7, our trapping mechanism provides a natural explanation for this initial 
conflguration. 

Horne and Moore have argued that the classical motion on certain moduli spaces is 



ergodic, provided that the potential energy is negligible. This means that all conflgurations 
are sampled given a sufliciently long time, and in particular a given modulus will eventually 
approach an ESP. We will argue that quantum corrections to the classical trajectory are 
signiflcant, and indeed lead to trapping, whenever the classical trajectory comes close to 
an ESP. Combining these two observations, we expect that in the full, quantum- corrected 
system the moduli are stuck near an ESP at late times. This means that the quantum- 
corrected evolution is not fully ergodic: the dynamics of (see also ||2^) implies that 



the modulus will eventually approach an ESP, at which point quantum effects will trap it 
there, preventing the system from sampling any further regions of moduli space. 



2. Moduli Trapping: Basic Mechanism 

We will now describe the mechanism of moduli trapping in more detail. Our discussion 
in this section will be based on simple estimates of particle production and the consequent 
backreaction, generalizing the results of to the case of a complex field. A more 

complete analysis, along with numerical results, will be presented in §3. 

We will consider the specific model 

L = ld^<Pd^4> + \d^xd^X - y l</'l V (2.1) 

where a complex modulus 4> = 4>i -\- i(f)2 interacts with a real scalar field x- We are 
restricting ourselves to the case of a fiat moduli space which has a single ESP at (/> = 0, 
where x becomes massless, and a particularly simple form for the x interaction. This 
simple case illustrates the basic physics and can be generalized as necessary, for example 
to include supersymmetry. 

We will consider the case where cp approaches the origin with some impact parameter 
following a classical trajectory of the form 

(f){t) =in + vt. (2.2) 

Classically, if x vanishes then ( |2.2| ) is an exact solution to the equations of motion, and 
the presence of the ESP will not affect the motion of 4>. 



Quantum effects will alter this picture considerably, because the trajectory (|2.2| ) will 
lead to the production of x particles, as we discuss in §2.1. The backreaction of these 
particles on the motion of (p will then lead to trapping, as we will see in §2.2. In §2.3 we 
will illustrate this effect with the example of colliding D-branes. 

2.1. Quantum Production of x Particles 

Let us first study the creation of x particles without considering how they may back- 
react to alter the motion of (p. In this approximation we may substitute ( |2.2[ ) into the 
action ( p.l|) to get a free quantum field theory for x with a time-varying mass 



ml{t)^g'\m'- (2.3) 
This time dependence leads to particle production. 
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Consider a mode of the x held with spatial momentum k, whose frequency 



u{t) = ^fWT^f^W)? (2.4) 

varies in time. This mode becomes excited when the non-adiabaticity parameter Cj/uj'^ 
becomes at least of order one. This parameter vanishes as t ^ ±00, indicating that 
particle creation takes place only while (f) is near the ESP. It is straightforward to see that, 
for the trajectory ( p.2|) , uj/uj'^ can be large only in the small interval |(^| < near the 
ESP, where 

A</> = (2.5) 

and only for momenta 

< 1. (2.6) 

gv 

When the quantity on the left hand side is small, particle creation effects are very strong. 
They are strongest if the modulus passes sufficiently close to the ESP, i.e. if 



In this case x modes whose momenta k fall in the range ( |2.6| ) will be excited.ll Qualitatively, 
we expect that the occupation numbers Uk of such modes will vary from zero (no real 
particles) for modes with vanishing non-adiabaticity to of order unity for modes with very 
large non-adiabaticity. The full computation of given in Appendix A yields 

n.=expf-7r^i±^y (2.9) 



which agrees with this qualitative expectation. Note that even when ( |2.7|) is not satisfied, 
there is generically a nonvanishing, though exponentially suppressed, number density of 
created particles; even in this case we will find a nontrivial trapping effect. 



^ This may be checked as follows. We have argued that unsuppressed particle production 
occurs only when the modulus is sufficiently close to the ESP, \(f>\ < ^Jvjg. The modulus remains 
within this window for a time 

At ~ ^^^^ - (5t;)"'/^ (2.8) 

The uncertainty principle implies in this case that the created particles will have typical energy 
E ^ (At)^^ and thus momenta k ^ {gv — g^fi^y^'^. This agrees with the estimate ( ^ 
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Before discussing the backreaction due to the production of x particles, it is crucial to 
control other effects from the x field. In particular, there is another important quantum 
effect which arises in motion toward the origin: loops of light x particles give corrections 
to the effective action. These include both kinetic corrections and the Coleman- Weinberg 
potential energy. The latter we will subtract by hand, as we will explain in §3.1. This gives 
a good approximation to the dynamics in any situation where kinetic energy dominates. 

The kinetic corrections are organized in an expansion in v'^ /(f)'^ P). The parameters 
controlling both remaining effects - the nonadiabaticity controlling particle production 
and the kinetic factor f^/^^ controlling light virtual x particles - diverge as we approach 
the origin. However, at weak coupling, the nonadiabaticity parameter is parametrically 
enhanced relative to the kinetic corrections, i.e. v"^ /g'^cj)^ ^ f^/^^, so we can sensibly 
focus on the effects of particle production. More specifically, we can ensure that the 
kinetic corrections are insignificant by including a sufficiently large impact parameter /x. 

We will also analyze the case of small /U, including fi = 0. This relies on the plausible 
assumption that the effects of the kinetic corrections remain subdominant as we approach 
very close to the origin, and that in particular in our weak coupling case they do not by 
themselves stop (p from progressing through the origin. It would be interesting to develop 
theoretical tools to analyze this issue more directly and check this hypothesis. 

2.2. Backreaction on the Motion of (p 

One might expect a priori that any description of the motion of (p which fully incorpo- 
rates backreaction from particle production would be immensely complicated. Fortunately, 
this turns out not to be the case, and a simple description is possible. The key simplifica- 
tion is that creation of x particles happens primarily in a small vicinity of the ESP ^ = 0, 
so one can treat this as an instant event of particle production. These particles induce 
a very simple linear, confining potential acting on (p, V ^ \(p\. The motion of (p in this 
potential between successive events of particle production can be described rather simply. 

Let us now explore this in more detail. We have seen that as (p moves in moduli 
space, some of its energy will be transferred into excitations of x- This leads to a quantum 
vacuum expectation value (x^) 7^ 0. As (p rolls away from the ESP, the mass of the created 
X particles increases, further increasing the energy contained in the x sector. At this point 
the backreaction of the x field on the dynamics of (p becomes important, and the moduli 
space approximation breaks down. 
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We will concentrate on the backreaction of the created particles on the motion of the 
field (j) far away from the small region of non-adiabaticity, i.e. for (f) ^ ~ g. At this 
stage the typical momenta are such that the x particles are nonrelativistic, k < ^Jgv <^ 
Therefore the total energy density of the gas of x particles is easily seen to be 

pM = I 0^n,^k^+g^\m'^9m)K, (2.10) 
where is the number density of x particles, 

As (p continues to move away from the ESP (p = 0, the number density of x particles 
remains constant, as particles are produced only in the vicinity of ^ = 0. However, the 
energy density of the x particles grows as g\(f){t)\n^. This leads to an attractive force of 
magnitude gn^, which always points towards the ESP (p = 0. 

This force of attraction slows down the motion of (p, and eventually turns (p back 
toward the ESP. This reversal occurs in the vicinity of the point (p^ at which the initial 
kinetic energy density = ^ matches the energy density contained in x particles. 
We find 

(P^ = i!!l^i/2e-9^V-. (2.12) 

Observe that for <^ 1 the trapping length on the first pass is always much greater than 
the impact parameter /x, which means that the motion of the moduli after the first impact 
is effectively one- dimensional. 

After changing direction at (p^, (p falls back toward the origin. On this second pass by 
the ESP, more x particles are produced, leading to a stronger attractive force. This process 
repeats itself, leading ultimately to a trapped orbit of (p about the ESP, in a trajectory 
determined by the effective potential and consistent with angular momentum conservation 
on moduli space. 

We conclude that, in this simplified setup, a scalar field which rolls past an ESP will 
oscillate about the ESP with an initial amplitude given by ( |2.12| ). 

In fact, in many cases the amplitude of these oscillations will rapidly decrease due to 
the effect of parametric resonance, similar to the effects studied in the theory of preheating 
1^ , and the field (p will fall swiftly towards the ESP. This important result will be described 
in more detail in 53.3. 
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So far we have not incorporated the effects of scattering and decay of the x particles. 
These could weaken the trapping potential (|2.10|) by reducing the number of x particles. 
Specifically, the energy density contained in a fixed number of x particles ( pj.lUD grows 
at late times, since the x mass increases as 4) rolls away from the ESP. However, if the 
number density of x particles decreases due to annihilation or decay into lighter modes, 
this mass amplification effect is lost. It is therefore important to determine the rate of 
decay and annihilation of the x particles. 

In Appendix B we address these issues and demonstrate that the trapping effect is ro- 
bust for certain parameter ranges, provided that the light states are relatively stable. This 
stability can easily be arranged in supersymmetric models, and in fact occurs automatically 
in certain D-brane systems. 

Rescattering effects, in contrast, may actually strengthen the trapping effect. Once x 
particles have been created, they will scatter off of the homogeneous 4> condensate, causing 



it to gradually decay into inhomogeneous, decoherent cj) excitations However, we 

will not consider this potentially beneficial effect here. 



2.3. The Example of Moving D-hranes 

Before proceeding, it may be illustrative to discuss these results in terms of a simple, 
mechanical example - a moving pair of D-branes. The moduli space of a system of two 
D-branes is the space of brane positions. In terms of the brane worldvolume fields the 
separation between the two branes can be regarded as a Higgs field (f). The off-diagonal 
components of the U (2) gauge field are the W bosons. At the ESP of this system, = 0, 
the W bosons are massless. Away from = the W bosons acquire a mass by the Higgs 
mechanism, breaking the symmetry group from U{2) down to U{1) x ^7(1). If we identify 
X with the W fieldi and g"^ ~ 9ym ~ 9s with the string coupling, then we find that the 
brane worldvolume theory contains a term like ( |2.1| ). We therefore expect this system to 
exhibit moduli trapping. 

^ For simplicity we ignore the superpartner of the x boson. 
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Fig. 1: This figure illustrates the creation of open strings as two D-branes pass 
near each other. The left corner shows the target space picture of the creation of 
the open strings. 

The trapping effect is a quantum correction to the motion of D-branes. As the D- 
branes approach each other, the open strings stretched between them become excited. 
When the D-branes pass by each other and begin moving apart the stretched open strings 
become massive and puU the D-branes back together. We depict this in fig. |1|. 

This effect can be a significant correction to the dynamics of any system with a 
number of mobile, mutuaUy BPS D-branes. Consider, for example, N D3-branes which fill 
spacetime and are transverse to a compact six-manifold M. Let us take these branes to 
begin with small, random, classical velocities in M. The classical dynamics of this system is 
similar to that of a nonrelativistic, noninteracting, classical gas. When we include quantum 
production of light strings, the branes begin to trap each other, pairwise or in small groups, 
then gradually agglomerate until only a few massive clumps of many branes remain. 

One interesting consequence is that such a system will tend to exhibit enhanced gauge 
symmetry, with gauge group U (N) if the final state consists of a single clump. (Hubble 
friction may bring the branes to rest before the aggregation is complete, in which case 
the gauge group will be a product of smaller factors; we will address related issues in 
§4.1.) Another important effect of massive clumps is their gravitational backreaction: a 
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large cluster of D-branes will produce a warped throat region in M, which may be of 
phenomenological interest [ p8| . 

There are additional corrections to the classical moduli space approximation of the 
D-brane motion which come from velocity-dependent forces. These correspond in the 
D-brane worldvolume field theory to higher-derivative corrections generated by virtual 
effects. When this field theory is at weak 't Hooft coupling, open string production is the 
dominant effect as one approaches an ESP. However, sufficiently large clusters of branes 
will be described by gauge theories at strong 't Hooft coupling, where the dynamics of 
additional probe branes is governed instead by the analysis of 

A similar interaction was studied in the context of the scattering of DO-branes in 



There is a crucial difference between that system and the case of interest here, in which the 
branes are extended along 34-1 dimensions. In the DO-brane problem, there is a nontrivial 
probability for the DO-branes to pass by each other without getting trapped: because the 
DO-brane is pointlike, there is some probability for no open strings between them to be 
created or for those created to annihilate rapidly. This is the leading contribution to the 
S-matrix. In our case, there is always a nonzero number density of particles created. As 
we argue in Appendix B, for certain ranges of parameters these particles do not annihilate 
rapidly enough to prevent trapping. 



3. Moduli Trapping: Detailed Analysis 

In the previous section we gave an intuitive explanation of the trapping effect, which 
we will now describe in more detail. In §3.1 we will present the equations of motion which 
govern the trajectory of the modulus (p, including the backreaction due to production of 
light particles. These equations are difficult to solve exactly, so in §3.2 we will integrate 
the system numerically. In §3.3 we focus on the special case /i = 0, where the modulus 
rolls directly through the ESP. In this case analytic techniques are available, and as we 
will see the trapping effect is considerably stronger than in the /U 7^ case. Even for /U 7^ 
some aspects can be derived analytically in a perturbative expansion. This last case is 
somewhat technical, so we refer the reader to Appendix C for details. 

^ A further correction to our dynamics could arise if, as we will discuss in §5, the branes keep 
moving until the system is beyond the range of effective field theory. 
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3.1. Formal Description of Particle Production Near an ESP 

The full equations of motion are found by coupling the classical motion of 4) to the 
time-dependent x quantum field theory defined by (|2.1|).i 

In general, the presence of an ESP will alter the moduli dynamics in two ways. First, 
any x excitations produced by the mechanism described above will backreact on the classi- 



cal evolution of (j). In particular, as we saw in ( |2.10|) , a non-zero expectation value (x^) 7^ 



arising from particle production effectively acts like a linear potential for cf) and drives the 
moduli towards the origin. This is the effect we wish to describe. Second, virtual x parti- 
cles generate quadratic and higher-derivative contributions to the effective action as well 
as an effective potential for a spacetime-homogeneous (f). 

As we discussed in §2.1, we can neglect the kinetic corrections in our weakly-coupled 
situation. The interaction in ( |2.1| ) also induces important radiative corrections to the 
effective potential. Specifically, it leads to a Coleman- Weinberg effective potential and 
three UV-divergent terms: 

V,ff{ct^) = A,ff+g'mlff<p'+g'X,ff<p\ (3.1) 

These UV divergences could be subtracted by hand using appropriate counterterms. In a 
supersymmetric system these divergences are absent. 

In order to isolate the effects of particle production at the order we are working, we will 
subtract by hand the entire Coleman- Weinberg effective potential for (p that is generated 
by one loop of x particles. This mimics the effect of including extended supersymmetry, 
which is a toy case of interest in string theory and supergravity. For the more realistic 
A/" = 1 supersymmetry in four dimensions, radiative corrections do generically generate 
a nontrivial potential energy. Nevertheless, particle production effects can still dominate 
the virtual corrections to the potential after spontaneous supersymmetry breaking. The 
reason is that bosons and fermions contribute with opposite signs in loops, but on-shell 
bosons and fermions, such as those produced by the changing mass of x, contribute with 
the same sign to backreaction on (p. 

To describe the production of x particles, we first expand the quantum field x in terms 
of Fock space operators as 

X = ^akXk + alxl (3.2) 



We remain in flat space quantum field theory, reserving gravitational eflFects for §4. 
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where the Xk ai'e a complete set of positive-frequency solutions to the Klein-Gordon equa- 
tion with mass 

ml{t)^g'\m'- (3.3) 



Expanding in plane waves 
the equation of motion is 



Xfc = txfc(t)e'"-" (3.4) 



(a2 + A;2 + /|^(t)|2jwfc = 0. (3.5) 

The modes ( p.4|) are normalized with respect to the Klein-Gordon inner product, which 
fixes 

uliik - iiluk = -i. (3.6) 



The wave equation ( |3.5| ) has two linearly-independent solutions for each /c, so in general 
there will be many inequivalent choices of positive- frequency modes Xk- Each such choice 
of mode decomposition defines a set of Fock space operators via ( p.2|) , which in turn define 
a vacuum state of the theory. The wave equation depends explicitly on time, so there 
is no canonical choice of Poincare invariant vacuum. Instead, there is a large family of 
inequivalent vacua for x- 

We can choose a set of positive frequency modes w^" that take a particularly simple 
form in the far past. 



-^r^k^+,^^it'Wdt' ^,t^_^_ (3.7) 



This choice of mode decomposition defines a vacuum state \in). In the far past the phases 
of the solutions ( p.7|) are monotone decreasing with t, indicating that the state \in) has no 
particles in the far past. This state, known as the adiabatic vacuum, evolves into a highly 
excited state as the modulus (j) rolls past the ESP. 

We can now write down the classical equation of motion for (f) including the effects of 
X production. Including a subtraction 5m , to be determined shortly, it is 

{d^ + g\{x^)-5M)y = ^. (3.8) 

The expectation value (x^) depends on time and is calculated in the adiabatic vacuum 
\in). At time t 

{tn\x\t)\in)= / j—^\ur{t)\\ (3.9) 



(27r) 
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where the u^^ are determined by the boundary condition (|3.7| ) in the far past. 

In order to subtract the Coleman- Weinberg potential, we must remove the contribution 
to (x^) coming from one loop of x particles, replacing the x mass-squared with g'^\(p(t)\'^ . 
That is, the subtraction Sm can be written as 

5m = f ,t\ , ^ - (3.10) 

} (27r)3 2^/FTFH' 

With this form it is straightforward to see that when the impact parameter is very large, 

((x^) ~ ^m) is negligible and ^ follows its original trajectory (|2.2| ). 

To summarize, the effects of quantum production of x particles on the classical motion 

of the modulus ^ are governed by: 



{di + e+g'\m?yi:-^ (3.11) 



^^'^'^^ = /(0i"™i'- 

The above equations of motion can be reformulated in terms of the energy transferred 
between the two systems. In particular, it is straightforward to show that the coupled 
equations (|3.11| ) are equivalent to the statement 

|^<^ = -^(^^l^xN^)- (3-12) 

The left-hand side of (|3.12|) involves the classical energy of the rolling (^(t) fields, whereas 
the right hand side is an expectation value of the time-dependent x Hamiltonian calculated 
in quantum field theory. This is the more precise form of energy conservation which applies 
to our rough estimate in §2.2. 

Furthermore, the angular momentum on moduli space is conserved, since the action 
( p.l|) is invariant under phase rotations 4> ~^ 4>e^^ . In the present case (|2.1|) , the x particles 
do not carry angular momentum, so the orbit of (j) around the ESP will have fixed angular 
momentum. The result is an angular momentum barrier which keeps the modulus at a 
finite distance from the ESP. 

More complicated scenarios allow for the exchange of angular momentum between (p 
and X- This includes the case of colliding D-branes, where the strings stretching between 
the two D-branes can carry angular momentum. Moreover, as we will see in §4, the situ- 
ation changes once gravitational effects are included, as angular momentum is redshifted 
away by cosmological expansion. This leads to scenarios where the moduli are trapped 
exactly at the ESP, rather than orbiting around it at some finite distance. 
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3.2. Moduli Trapping: Numerical Results 

The coupled set of integral and differential equations (|3.11|) governing the trapping 
trajectory is hard to solve in general. Some analytic results can be obtained through an 
expansion in the non-adiabaticity parameter ujjup' ., combined with a systematic iteration 
procedure. Specifically, to the extent that this non-adiabaticity parameter is small enough 
that the correction to the motion of the moduli field only shows up at exponentially late 
times, we can calculate the bending angle and the energy loss during its first pass. This is 
given in Appendix C. 

As time goes on, the mass amplification of the x particles makes higher order terms 
as well as non-perturbative terms in the adiabatic expansion crucial for the motion of the 
moduli. This makes it very hard to proceed analytically to obtain the detailed evolution 
of the system. 

We have numerically integrated the coupled equations ( |3.11|) in Mathematica, using a 
discrete sum to approximate the momentum integral /c, and implementing the subtraction 
of the Coleman- Weinberg potential described above. 




Fig. 2: This figure shows the evolution, in the complex (\) plane, of a system with 
parameters = 20, /x = 0.3, = 1. The field rolls in from the right and gets 
trapped into the precessing orbit exhibited in the plot. The orbit is initially an 
elongated ellipse, but gradually becomes more circular. In an expanding universe, 
the field would lose its angular momentum, so that the radius of the circle would 
eventually shrink to zero. 

In fig. 1^ we plot a trajectory for the case /i > 0, where (j) becomes trapped in a 
spiral orbit around the ESP. The radius of the orbit varies with the parameters, but the 
qualitative features shown are typical. 
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Fig. 3: This shows one-dimensional trapping, in which </> passes directly through 
the ESP (f) = 0. The vertical axis is the real part of (j), and the horizontal axis 
is time. The amplitude of the oscillations decreases exponentially as a result of 
parametric resonance, as we explain in §3.3. 

In fig. ^ we plot the trajectory of a modulus which is aimed to pass directly through 
an ESP, with vanishing impact parameter. In this case the motion becomes effectively 
one-dimensional, and the field moves directly through the ESP ^ = 0. The trapping effect 
in this case is especially strong, and can be understood analytically to come from resonant 
production of x particles, as we will now explain. 

3.3. The Special Case of One- Dimensional Motion 

In this section we will concentrate on the interesting and important special case of 
one-dimensional motion, i.e. vanishing impact parameter fi. Perhaps surprisingly, this 
is a good approximation to the general case. Indeed, the results of §2 demonstrate that 
trapping becomes exponentially suppressed when the impact parameter (the imaginary 
part of the moduli field) becomes greater than On the other hand, for /i <^ 



the motion of the field (p stops at 4>^ ~ ^5^3 ■ The ratio of to /x in the regime where 
trapping is efficient (i.e. for n < \/ ^) is therefore 



Thus, in the case of efficient trapping and weak coupling, the ellipticity of the moduli orbit 
is very high, so that the motion is effectively one-dimensional. 







(3.13) 
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In the case fi = the number density of x particles created when the field cp passes 
the ESP is 



n = (3.14) 

(27r)3 ^ ' 

At 1^1 ^ ^J^, when the x particles are nonrelativistic, the mass of each particle is equal 
to g\(j)\, and their energy density is given by 

pM) = = \<f>\ (3.15) 

We have written \(p\ because this energy does not depend on the sign of the field (p. This 
will be very important for us in what follows. 

One should note that, strictly speaking, the x particles have some kinetic energy even 
at (p = 0, but for fi* <^ 1 this energy is much smaller than the kinetic energy of ^ [0]: 

pJ'^=«)~S^."^-^pT- (3.16) 

This means that the energy of (p decreases only slightly when it passes through the ESP 
(p = 0. Although the initial energy in x particles is small, this energy increases with \(p\, 
~ gn^\cp\, and creates an effective potential for cp. The equation of motion for cp in this 
potential is [|1|: 

-<P + gn^^^=0. (3.17) 

The last term means that is attracted to the ESP ^ = with a constant force proportional 
to n^. 

At some location cp\ the x energy density equals the initial kinetic energy density 
= it)^; at this point (p stops and then falls back toward ^ = 0. 

On this second pass by the origin, the energy density of the x particles again becomes 
much smaller than the kinetic energy of (p. Energy conservation implies that (p will pass 
the point (/> = at almost exactly the initial velocity v. Since the conditions are almost the 
same as on the first pass, new x particles will be created, i.e. will increase. The field 
will continue moving for a while, stop at some point cp2-, and then fall back once more 
to the ESP, creating more particles. Because each new collection of particles is created 
in the presence of previous generations of particles, the process occurs in the regime of 
parametric resonance, as in the theory of preheating. 

18 



A detailed theory of this process was considered in [^]; see in particular Eqs. (59), (60). 
By translating the problem into a one- dimensional quantum mechanics system (as in Ap- 
pendix A) with a particle scattering repeatedly across an inverted harmonic potential, 
1^ calculated the multiplicative increase of the Bogoliubov coefficients during each pass 
in terms of the reflection and transition amplitudes. In application to our problem, the 
equations describing the occupation numbers of x particles with momentum k produced 
when the field passes through the ESP j -f 1 times look as follows: 

nf^=n^,exp(2Vfc), (3-18) 

where 

^^fc = ^ln(^l + 2e-^^' -2sin^-'' e"^^' V^TTe^^) . (3.19) 
Here i'^ = ^ and 6^ is a relative phase variable which takes values from to 2n. In a 

gv ^ 

cyclic particle creation process in which the parameters of the system change considerably 
during each oscillation (which is our case, as will become clear shortly), the phases 6^ 
change almost randomly. As a result, the coefficient iij for small k takes different values, 
from 0.28 to —0.28, but for 3/4 of all values of the angle 6^ the coefficient iij is positive. 
The average value of /ij is approximately equal to 0.15. This means that, on average, the 
number density of x particles grows by approximately a factor of two or three each time 
that (j) passes through the ESP = 0. 

But this means that with each pass, the coefficient in ( |3.15| ) grows by a factor 
of two or three. It follows that the effective potential becomes two to three times more 
steep with each pass. Correspondingly, the maximal deviation |^*| from the point (/> = 
exponentially decreases with each new oscillation. Since the velocity of the field at the 
point = remains almost unchanged until cj) loses its energy to the created particles, the 
duration of each oscillation decreases exponentially as well. Therefore the whole process 
takes a time 0{lQ)(f)\/v, after which the backreaction of the created particles becomes 
important, and the field falls to the ESP. 

This process is very similar to the last stages of preheating, as studied in [|1|. The 
main difference is that in the simplest models of preheating the field oscillates near the 
minimum of its classical potential. In our case the effective potential is initially absent, 
but a potential is generated due to the created particles. This is exactly what happens at 
the late stages of preheating, when the effective potential (with an account taken of the 
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produced particles) becomes dominated by the rapidly-growing term proportional to |(^|; 
see the discussion in Section VIII B of [0. 

We would like to emphasize that until the very last stages of the process, the backre- 
action of the created particles can be studied by the simple methods described above. At 
this stage the total number of created particles is still very small, but their number grows 
exponentially with each new oscillation. This leads to an exponentially rapid increase of 
the steepness of the potential energy of the field (j) ( |3.15|) and, correspondingly, to an expo- 
nentially rapid decrease of the amplitude of its oscillations. This extremely fast trapping 
of (j) happens despite the fact that at this first stage of oscillations the total energy of (j), 
including its potential energy, remains almost constant. 

Once the amplitude of oscillations becomes smaller than the width of the nonadi- 



abaticity region, < ~ one can no longer assume that the number of 

particles will continue to grow via a rapidly-developing parametric resonance. The ampli- 
tude of the oscillations is given by 2g\4>\n ' amplitude becomes 0{\Jv/ g) when the 
total number of the produced particles grows to 



~ ^3/2^-1/2 _ (3 20) 



Note that the typical energy of each x particle at ^ yv/g is oi the same order as 

its kinetic energy 0{^/gv). One can easily see that the total energy density of particles x 
at that stage is roughly yfgvn^ ~ i.e. it is comparable to the initial kinetic energy 

of (j). 

Thus, our estimates indicate that the regime of the broad parametric resonance ends 
when a substantial part of the initial kinetic energy of (f) is converted to the energy of the 
X particles, and the amplitude of the oscillating field (j) becomes comparable to the width 
of the nonadiabaticity region, 

|(/>| ~ A(/) = . (3.21) 

We will use these estimates in our discussion of the cosmological consequences of 
moduli trapping. In order to obtain a more complete and reliable description of the 
last stages of this process one should use lattice simulations, taking into account the 
rescattering of created particles [ p6| , |27[1 . An investigation of a similar situation in the 



theory of preheating has shown that rescattering makes the process of particle production 
more efficient. This speeds up the last stages of particle production and leads to a rapid 
decay of the field (j) HHl ? which in our case corresponds to a rapid descent of (f) toward the 
enhanced symmetry point. 
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4. Trapped Moduli in an Expanding Universe 



4-1- Rapid Trapping 

In this section we will study tiie conditions under wliicli tlie trapping meclianism in 
quantum field theory survives the effects of coupling to gravity in an expanding universe. 

First, we should point out one very beneficial effect of cosmological expansion. The 
field-theoretic mechanism presented above often leads to moduli being trapped in large- 
amplitude fluctuations ( |2.12| ) around an ESP when fi ^ 0. On timescales where the 
expansion is noticeable, Hubble friction will naturally extract the energy from this motion, 
drawing the modulus inward and leading the modulus to come to rest at the ESP. 

Let us now ask whether the expansion of the universe can impede moduli trapping. 
Consider a system of moduli coupled to gravity, with the fields arranged to roll near an 
ESP. For simplicity we will consider FRW solutions with flat spatial slices, 

ds'^ = dt^ - a{tfdx^. (4.1) 



The Friedman equation determining a{t) is 

= (4.2) 

p 

where H = d/a and p is the energy density of the moduli. 

The trapping effect will be robust against cosmological expansion if the timescale 
governing trapping is short compared to H~^, i.e. if H <^ f where (p^ is given by 
( p.l2|) . Assuming that the potential energy of the moduli is non-negative, this implies that 

« ^TTi-^^e-- « 1. (4.3) 

This condition suffices to ensure that trapping is very rapid. 

If this condition is satisfied, trapping occurs in much less than a Hubble time, in which 
case the analysis of §2 and §3 remains valid. We will show in §4.3 that even when ( [4.3| ) is 
not satisfied, trapping does still occur, although with somewhat different dynamics. 
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4-2. Scanning Range in an Expanding Universe 

An important effect of the gravitational coupling is that during the expansion of the 
universe, the energy density in produced x particles dilutes like if they are non- 

relativistic and like if they are relativistic. The energy in coherent motion of (p, 

however, has the equation of state p = p and therefore dilutes much faster, as 

This effect reduces the range of motion for the moduli even before they encounter any 
ESPs. Hubble friction slows the progress of any rolling scalar field, and if the distance 
between ESPs is sufficiently large then a typical rolling modulus will come to rest without 
ever passing near an ESP. In order to apply our results to the vacuum selection problem, 
we will need to know how large a range of (p we can scan over in the presence of Hubble 



friction. This can be obtained as follows ||3C 
If we are in an FRW phase, 

a{t) = aot^ (4.4) 
then the equation of motion for (p (ignoring any potential terms) 

^ + 3H(i) = (4.5) 



has solutions of the form 



to 



<j>{t) = v{^] . (4.6) 



We can integrate this to determine how far the field rolls before stopping. 

Let us first consider the case (3 = 1/3, which corresponds to the equation of state 
p = p. This includes the case where the coherent, classical kinetic energy of (p drives the 
expansion. The value of (p, 

4>{t) = vto log (f) (4.7) 



diverges at large t. Thus (p can travel an arbitrarily large distance in moduli space. 
In the more general case /3 > 1/3 the field will travel a distance 

before stopping. 

In order to be in a phase with f3 > 1/3, the kinetic energy of (p must not be totally 
dominant; that is, we must have |(/)^ < p, where p = 3MpH^ is the total energy density 
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appearing on the right hand side of the Friedman equation. Plugging this into ( [4 .81) we 
obtain the constraint 

cp(t) - (Pito) < V6Mp^^. (4.9) 

Let us consider a specific example. Suppose that we start at to with kinetic energy 
domination: Kq/pq = 1 — e, e ^ 1, in some region of the universe that can be modelled as 
an expanding FRW cosmology. The kinetic energy drops like K ~ po(oo/a)^ ~ po(^o/^)^, 
while the other components of the energy dilute like 

p(t)=e/9o(toA)'+", (4.10) 

with w < 1. The universe will stop being kinetic-energy dominated at the time tc — 
tQe~^/'^^~'^\ at which point, according to ([4.7|), the modulus has travelled a distance 

(l){t^)-(l){to) = -—^vtQ\oge. (4.11) 
1 — w 

After this the field keeps moving and covers an additional range 

(t^{U) - <Pitc) = \/3M / (4.12) 

6{1 — W) 

To get a feel for the numbers, consider the case where vto ~ Mp, e ~ 10~^, and 
w = 0. Then cp will travel a total distance (p{t^) — (p{to) ~ 6Mp in field space, which is not 
particularly far. However, as we will discuss in §6, certain moduli spaces of interest have 
a rich structure on sub-Planckian scales, so in these cases there is a good chance that the 
modulus will encounter an ESP and get trapped before Hubble friction brings the system 
to rest. 

There is another natural possibility if we assume low-energy A/" = 1 supersymmetry. 
If the moduli acquire their potentials from supersymmetry breaking then there is a large 
ratio between the Planck scale and the scale of these potentials, leading to significant 
scanning ranges. Specifically, consider a contribution to the energy density coming from 
a potential energy V at the supersymmetry-breaking scale. If the initial kinetic energy of 
the moduli is Planckian and the supersymmetry-breaking scale is TeV then there will be a 
prolonged phase in which kinetic energy dominates, since e = V/M^ ~ 10~^^. This allows 
(p to scan a significantly super-Planckian range in field space. 
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4-3. Trapping in an Expanding Universe 

We are now in a position to combine all the relevant effects and consider trapping 
during expansion of the universe. For simplicity, we will concentrate on the case of ef- 



fectively one-dimensional motion, ^ <^ y^/g. Suppose that, taking into account Hubble 
friction, the modulus field passes in the vicinity of the ESP at some moment to, so that 
X particles are produced, with n-)^(to) = ^^^^^3 . We will now determine the remaining 
evolution including both our trapping force and Hubble friction. After the particles have 
been produced, the field becomes attracted toward = by a force gn^, so taking into 
account the dilution of the produced particles, for > the equation of motion is 

4> + 3H^ = -gn^{to)(^^y (4.13) 
For the general power law case, a{t) oc t^, this becomes 

+ 3^0 = -gn^ito)ito/tf^ • (4.14) 
The general solution of this equation is 

0(t) = 0(to) + c{t-^'^' - t-'^+') + - 0§^ito/tr't- (4.15) 

where c is some constant. In the important case (3 — 2/3, which corresponds to a universe 
dominated by pressureless cold matter, the general solution is 

<l) = m + c{t^' -t~')-9nxtllog^ . (4.16) 

to 

According to these solutions, in a universe dominated by matter with non-negative pressure 
(i.e. (3 < 2/3) the field (/) moves to —00 as t — > cxd. 

Of course, as soon as the field reaches the point (p = 0, this solution is no longer 
applicable, since the attractive force changes its sign (the potential is proportional to |0|). 
The result above simply means that the attractive force is always strong enough to bring 
the field back to the point (p — within finite time. Then the field moves further, with 
ever decreasing speed, turns back again, and returns to = once again. The amplitude 
of each oscillation rapidly decreases due to the combined effect of the Hubble friction and 
of the (weak) parametric resonance. This means that once passes near the ESP, its fate 
is sealed: eventually it will be trapped there. 



24 



4-4 ■ Efficiency of Trapping 

It is useful to determine what fraction of all initial conditions for the moving moduli 
lead to trapping. There are several constraints to be satisfied. First of all, if the impact 
parameter n is much larger than \/vJg, the number of produced particles will be expo- 
nentially small, and the efficiency of trapping will be exponentially suppressed. Of course, 
eventually (p will fall to the enhanced symmetry point, but if this process takes an expo- 
nentially large time, the trapping effect will be of no practical significance. Thus one can 



roughly estimate the range of interesting impact parameters to be Oi^^fvjg). 

Another constraint is related to the fact that even if initially the energy density of the 
universe was dominated by the moving moduli, as discussed in §4.2, these fields can only 
move the distance given by ( [4.11|) ,( ^.12 ). This distance depends on the initial ratio 1 — e 



of kinetic energy to total energy, leading to a scanning range CM^ in field space, where 
the prefactor C is logarithmically related to e. 

Thus, the field becomes trapped only if there is an enhanced symmetry point inside a 



rectangle with sides of length CMp along the direction of motion and width 0{^Jv/ g) in 
the direction perpendicular to the motion. 

Interestingly, the total area (phase space) of the moduli trap 

Strap -CMp^l (4.17) 

increases as the coupling decreases. This implies that the efficiency of trapping grows at 
weak coupling. Although this may seem paradoxical, it happens because the mass of the 
X particles is proportional to the coupling constant and (fixing the other parameters) it is 
easier to produce lighter particles. On the other hand, if g becomes too small, the trapping 
force gn^ ~ (^5/2^3/2 bg^omes smaller than the usual forces due to the effective potential, 
which we assumed subdominant in our investigation. 

So far we have studied the simplest model where only one scalar field becomes massless 
at the enhanced symmetry point. Let us suppose, however, that N fields become massless 
at the point ^ = 0. If these fields interact with (f) with the same coupling constant g, 
then particles of each of these fields are produced, and the trapping force becomes N times 
stronger. In other words, the trapping force is proportional to the degree of symmetry at 
the ESP. 
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5. String Theory Effects 

It is interesting to ask if there is any controUed situation where string-theoretic effects 
become important for moduh trapping. Here we wiU simply hst several circumstances 
in which stringy and/or quantum gravity effects can come into play, as well as some 
constraints on these effects, leaving a full analysis of this subtle and interesting situation 
for the future. 

5.1. Large x Mass 

One way stringy and quantum gravity effects could become important in the colliding 
D-brane case is if the x mass at the turnaround point is greater than string scale, g(p^ > mg. 
This can happen even if the velocity is so small that during the non-adiabatic period near 
the origin only unexcited stretched strings are created. Then, as in our above field theory 
analysis, we have 

9<P* = ^v'/\-^^'/\ (5.1) 

In this case, the full system includes modes, namely the created x strings, which are 
heavier than the string oscillator mode excitations on the individual branes. This means 
that the system as a whole cannot consistently be captured by pure effective field theory. 
However, it may still happen that the created stretched strings are relatively stable against 
annihilation or decay into the lighter stringy modes. Their annihilation cross section is 
suppressed by their large mass, as discussed in Appendix B.i Furthermore, an individual 
stretched string will not directly decay if it is the lightest particle carrying a conserved 
charge. 

This latter situation happens in the simplest version of a D-brane collision. The 
created stretched string cannot decay into lighter string or field theory modes because it 
is charged and they are not. 

5.2. Large v and the Hagedorn Density of States 

If we increase the field velocity (p = v, then we may obtain a situation in which excited 
string states are produced as (p passes the ESP. The number of string states produced in 

^ For stringy densities of stretched strings, there could be additional corrections to the anni- 
hilation rate, but we will not consider this possibility. 
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this process is enhanced by the Hagedorn density of states, so the Bogohubov coefficients 
have the structure 

n 

where in the D-brane context, g = ^Jgl is the Yang-Mills coupling on the D-branes. Because 
of the e~'^'^'^ I ^^""^ suppression in the second factor, this effect is only significant if gv ^ m^. 

However, in the case of colliding D-branes, and any situation dual to it, there is a 
fundamental bound on the field velocity from the relativistic speed limit of the branes. 
That is, for large velocity one must include the full Dirac-Born-Infeld Lagrangian for 
which takes the form 



^ = - fo '2 I 1 - (5-3) 

{2'n:ygsa'^ J \ m'l 

This action governs the nontrivial dynamics of (p for velocities approaching the string scale, 
and in particular, it refiects the fact that the brane velocity gcpa' must be less than the 
speed of light in the ambient space. Applied to our situation, ( ^.3|) implies that the D- 
brane velocity cannot be large enough for the Hagedorn enhancement ( ^.2|) to substantially 
increase the trapping effect. 

However, in the presence of a large velocity, the effective mass of the stretched string 
also has important velocity-dependent contributions This may enhance the non- 

adiabaticity near the origin and thus enhance the particle production effect. If anything, 
we expect this to increase the trapping effect; it would be interesting to study this case 
further. As we discuss at the end of Appendix A, by using unitarity combined with the 
stringy calculations in PH], it might be possible to determine the net result of all these 



effects on the stringy trapping mechanism. It would be interesting if, taking into account 
all relevant quantum effects, a large contribution to string production occurred in some 
controlled setting, as suggested in pl| , p2| , p3| . If such an effect occurred for motion near 



an ESP, it would apparently enhance the trapping effect and possibly even indicate that (p 
slows down enormously before passing through the origin, as occurred in the case studied 
in 0. However, it would be important to check for control of the quantum corrections to 
such a system. 
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5.3. Light Field- Theoretic Strings 

A further possibility is to formally reduce the tension of strings by considering strings 
in warped throats, strings from branes partially wrapped on shrinking cycles, and the 
like. In these situations, the strings are essentially field-theoretic, though string theory 
techniques such as AdS/CFT and "geometric engineering" of field theories may provide 
technical help in analyzing the situation. 



6. The Vacuum Selection Problem 

We can now apply the ideas of the previous five sections to the cosmology of theories 
with moduli. 

A natural application of the moduli trapping effect is to the problem of vacuum 
selection. One mechanism of vacuum selection is based on the dynamics of light scalars 
during infiation. Moduli fields experience large quantum fiuctuations during infiation and 
can easily jump from one minimum (or valley) of their effective potential to another. It 
was suggested long ago that such processes may be responsible, e.g., for the choice of 
the vacuum state in supersymmetric theories and for the smallness of the cosmological 
constant . The probability of such processes and the resulting field distribution depends 



on the details of the infiationary scenario and the structure of the effective potential . 

The mechanism that we consider in this paper is, in a certain sense, complementary 
to the infiationary mechanism discussed above. During infiation the average velocities of 
the fields are very small, but quantum fiuctuations tend to take the light scalar fields away 
from their equilibrium positions. On the other hand, after infiation, the fields often find 
themselves not necessarily near the minima of their potentials or in the valleys correspond- 
ing to the fiat directions, but on a hillside. As they roll down, they often acquire some 
speed along the valleys, see e.g. . At this stage (as well as in a possible pre- infiationary 
epoch) the moduli trapping mechanism may operate. 

This mechanism may reduce the question of how one vacuum configuration is selected 
dynamically out of the entire moduli space of vacua to the question of how one ESP is 
selected out of the set of all ESPs. This residual problem is much simpler because ESPs 
generically comprise a tiny subset of the moduli space. 
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6.1. Vacuum Selection in Quantum Field Theory 

In pure quantum field theory, discussed in §2, we saw that if a scalar field (p is initially 
aimed to pass near an ESP, then (p gets drawn toward the ESP and is ultimately trapped 
there. This appears to be a basic phenomenon in time-dependent quantum field theory: 
moduli which begin in a coherent classical motion typically become trapped at an ESP. 
This leads to a dynamical preference for ESPs. 

In many of the supersymmetric quantum field theories that have been studied rig- 
orously 0, the moduli space contains singular points at which light degrees of freedom 
emerge. We have seen that moduli can become trapped near these points given suitable 
initial conditions. 

6.2. Vacuum Selection in Supergravity and Superstring Cosmology 

Compactifications of M/string theory which have a description as a low energy effective 
supersymmetric field theory can have a natural separation of scales: the string or Planck 
scale can be much larger than the energy scales in the effective field theory potential. Thus, 
the intrinsically stringy effects of §5 are unimportant in this limit. On the other hand, the 
effects of coupling to gravity given in §4 continue to provide a crucial constraint, as we will 
now discuss. 

First of all, as in the case of pure quantum field theory, there exist very instructive 
toy models with extended supersymmetry, for which there is no potential at all on the 
moduli space. For these examples, in situations where higher-derivative corrections to 
the effective action are suppressed, a rolling scalar field has the equation of state p = p. 
This corresponds to the /3 = 1/3 case ( [4.7| ) of §4, for which one can scan an arbitrarily 
large distance in field space. Therefore, in this case, the trapping effect applies in a 
straightforward way to dynamically select the ESPs for regimes in which ( [4.3|) is also 
satisfied. 

More generally, however, one may wish to implement cosmological trapping in theories 
with some potential energy. In this case the requirement that the scanning range of (p (as 
constrained by Hubble friction in §4) should be large enough to cover multiple vacua is 
an important constraint. The absolute minimum requirement is that the scanning range 
is sufficient for the moduli to reach one ESP before stopping from Hubble friction; but to 
address the vacuum selection problem one should ideally scan a number of ESPs. 

One context in which this can happen is in a phase in which the kinetic energy of the 
rolling scalar fields dominates the energy density of the universe so that the (3=1/3 result 
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( [4 .71) applies. This may occur in a pre-inflationary phase in some patches of spacetime, 
though it is subject to the stringent hmitation in duration given in ( [4.11| ). Given such a 
phase, the field will roll around until it gets trapped at an ESP. 

During the ordinary radiation-dominated (/3 = 1/2) and matter-dominated (/3 = 2/3) 
eras, the more stringent constraint ( [4.9|) applies. As we indicated in §4, this scanning range 
is not large in Planck units, so we can usefully apply moduli trapping to the problem of 
vacuum selection in these eras only if the vacuum has appropriately rich structure on 
sub-Planckian scales. In other words, the average distance in moduli space between ESPs 
should be sub-Planckian. 

Gravitationally-coupled scalars (/> generically have a potential energy V{(l)/Mp) which 
has local minima separated by Planck-scale distances. In this cases, the limited scanning 
range during the /3 7^ 1/3 cosmological eras prevents our mechanism from addressing the 
vacuum selection problem. However, it is generic for compactification moduli to have 
special ESPs where the gravitationally-coupled system is enhanced to a system with light 
field theory degrees of freedom. Given a rich enough effective field theory in this ESP 
region, there will generically be interesting vacuum structure on sub-Planckian distances. 
In this sort of region moduli trapping will pick out the ESP vacua of the system. 

6. 3. Properties of the Resulting Vacua 

Let us now consider the qualitative features of the vacua selected by moduli trapping, 
assuming that the constraint imposed by Hubble friction has been evaded in one of the 
ways described above. 

First of all, it is important to recognize that what we have called ESPs may well be 
subspaces of various dimensions, not points. For example, in toroidal compactification of 
the heterotic string, there is one enhanced symmetry locus for each circle in the torus - new 
states appear when the circle is at the self-dual radius. Each of these loci is codimension 
one in the moduli space, but of course their intersections, where multiple radii are self-dual, 
have higher codimension. 

When moduli trapping acts in such a system of intersecting enhanced symmetry loci, 
we expect that the moduli will first become trapped on the locus of lowest codimension, 
but retain some velocity parallel to this locus. Further trapping events can then localize 
the modulus to subspaces of progressively higher codimension. The final result is that the 
moduli come to rest on a locus of maximally enhanced symmetry. 



30 



The simplest examples of this phenomenon are toroidal compactification, in which all 
circles end up at the self-dual radius, and the system of D-branes discussed in §2.3, in 
which the gauge symmetry is enhanced to [/(A^).0 

Quite generally, we expect that within the accessible range in field space, taking into 
account Hubble friction and the form of the potential, moduli trapping will select the ESPs 
with the largest number of light states, which often corresponds to the highest degree of 
symmetry.! 

In some very early epoch the rolling moduli can have large velocities, so trapping 
can occur even at points where the "light" states x have a relatively large mass, and the 
enhanced symmetry is strongly broken. However, Hubble friction inevitably slows the 
motion of the moduli. Thus, trapping at late times is possible only at ESPs with weakly- 
broken symmetries and very light particles. One could speculate about a possible relation 
of this fact to the mass hierarchy problem. 

Note that even though we emphasized the natural role of enhanced symmetry in 
moduli trapping, in fact the only strict requirement was the appearance of new light 
particles at the trapping points. In some of the many vacua of string theory, particles may 
be light not because of symmetry but because of some miraculous cancellations. Invoking 
such unexplained cancellations to produce a small mass is highly undesirable. However, 
moduli trapping may ameliorate this problem, as those rare points in moduli space where 
the cancellation does happen are actually dynamical attractors. 

Thus, the attractive power of symmetry and of light particles may have implications 
for questions involving the distribution of vacua in string theory P7|JH^ , P5| , |^ . Given the 



strong preference we have seen for highly-enhanced symmetry, the distribution of all string 
vacua obtained by a naive counting, weighted only by multiplicity, may be quite different 
from the distribution of vacua produced by the dynamical populating process discussed in 
our paper. It is therefore very tempting to speculate that some of the surprising properties 
of our world, which might seem to be due to pure chance or miraculous cancellations, in 
fact may result from dynamical evolution and natural selection. 



63P 



^ A toy model for this situation, in the case of three D-branes, has the potential ^ [xi|</'2 — 
+ xil'/'i — </'3p + xil'/'i ~ "J^ap], where (f)i and Xi are six different fields. Suppose that 02 
moves through the point (/)2 — (/'a = 0. This creates xi particles and traps the system at 02 = 03, 
where xi is massless. Subsequent motion of 0i can trap it at the point 0i = 02 = 03, making the 
remaining fields X2 and X3 massless. 

^ Moreover, as we discuss in Appendix B, the trapping effect is far more effective at ESPs for 
which the x particles do not decay rapidly. We therefore expect moduli trapping to select ESPs 
which have relatively stable light states. 
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7. The Moduli Problem 



One aspect of the moduU problem is that reheating and nucleosynthesis can be cor- 
rupted by energy locked in oscillations of the moduli. The source of the problem is that 
the true minima of the low-temperature effective potential applicable after inflation do 
not coincide with the minima of the Hubble-temperature effective potential which is valid 
during inflation. It follows that moduli which sit in minima of the latter during inflation 
will find themselves displaced from their true, low-temperature minima once inflation is 
complete. The energy stored in this displacement, and in the resulting oscillations about 
the true minimum, poses problems for nucleosynthesis. 

One way to address this problem is to permit initial displacements of the moduli, as 
described above, but somehow arrange that the oscillating moduli decay very rapidly to 
Standard Model particles. Alternatively, one could fix the moduli at a scale high enough 
that the Hubble temperature during inflation does not destabilize them. This may work 
in string models with stabilized moduli such as [^ , p| , p9t^ |. 



Another approach to this problem [|TT| is to posit that the moduli sit at an enhanced 
symmetry point minimum of the finite-temperature effective potential during inflation. 
Then, when inflation ends, the moduli are still guaranteed to be at an extremum of the 
effective potential. If this extremum is a minimum then the moduli have no problematic 
oscillations after inflation. Our trapping mechanism allies nicely with this idea by providing 
a preinflationary dynamical mechanism which explains the initial condition assumed in this 
scenario. That is, in parts of the universe where (j) kinetic energy dominates well before 
inflation, the trapping effect can explain why the moduli find themselves in ESP minima 
at the onset of infiation. 



8. Trapped Infiation and Acceleration of the Universe 

The main motivation of our investigation was to study the behavior of moduli in 
quantum field theory and string theory. However, the results we have obtained have more 
general applicability. To give an example, in this section we will study the cosmological 
implications of the trapping of a scalar field (f) with a relatively steep potential. 

Consider the theory of a real scalar field (j) with the effective potential m^(^^/2. In 
the regime cj) < Mp the curvature of the effective potential is greater than H^, with H 
the Hubble parameter, so (p falls rapidly to its minimum, and infiation does not normally 
occur. 
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We will assume that (j) gives some bosons x ^ mass g\(t) — (t)i\. Let us assume that 4> 
falls from its initial value cj)Q = -\- a) < Mp with vanishing initial speed. If we take 
a -C 1 and neglect for the moment the expansion of the universe, then (p arrives at 4>i with 
the velocity v = \pla mcpi . 

As (p passes (pi, it creates x particles with number density = (gv)^^"^ /8n^. After a 
very short time these particles become nonrelativistic, and further motion of (p away from 
(pi requires an energy g\(p — (pi\n^. In other words, the effective potential becomes 

V{(P) ^ lm'(P^+gn^\(P-(Pi\ = ]W(P^ + g'>/^'^\(P - (Pi\. (8.1) 

For g^/"^ STT^m^ ^ minimum of the effective potential is not at = 0, but at the 

point (pi, where the particle production takes place. The condition g^^'^-^^r^ > (pi implies 
that 

m < 2-^/\-^g^a^/^(Pi. (8.2) 

Thus, if the mass of (p is sufficiently small, the field will be trapped near the point (pi. 
To give a particular example, take (pi ~ Mp/2, a ~ 1/4. Then (p is trapped near (pi if 

m < 10"V^p- (8.3) 

For a very light field, such as a modulus with m ~ 10^ GeV ~ 10~^^Mp, this condition is 
readily satisfied unless g is very small. 

Once the field is trapped, it starts oscillating around (pi with ever-decreasing am- 
plitude, creating new x particles in the regime of parametric resonance. Eventually (p 
transfers a large fraction of its energy to x particles. One can easily check that in this 
model the fall of (p to the point (pi and the subsequent process of creation of x particles 
occurs within a time smaller than H~^, so one can neglect expansion of the universe at 
this stage. This process is therefore governed by the theory described in §3. In particular, 
we may use the estimate (|3.20| ) of the total number of x particles produced in the process. 
At the end of the particle production, the correction to the effective potential becomes 
much larger than at the beginning of the process: 

AV = g\(P - (Pi\n^ ~ v^/^g^/^\(P - (/>i| , (8.4) 

Subsequent expansion of the universe dilutes the density of x particles as a~^, which 
eventually makes the correction to the effective potential small, so that (p starts moving 
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down again. The field (f) remains trapped at (j) = (pi until the scale factor of the universe 
grows by a factor 



since the beginning of the trapping process. 

In the beginning of the first e-folding, the kinetic energy of the x particles and of 
the oscillations of (p is comparable to the potential energy of (p. However, the kinetic 
energy rapidly decreases, and during the remaining time the energy is dominated by the 
potential energy V{(f)i). This means that the trapping of (p may lead to a stage of inflation 
or acceleration of the universe, even if the original potential V{(p) is too steep to support 
inflation. 

Let us consider various possibilities for the scales in the potential, to get some simple 
numerical estimates for the duration of inflation. For example, if we take a,g = 0{1), 
(pi ~ Mp and m ~ 10^ GeV, then the scale factor during a single trapping event will 
grow by a factor of e^. If one considers a model with m ~ 10-3°Mp, which can arise m a 
radiatively stable manner (as in the "new old inflation" model fl^), the scale factor during 
a single trapping event can grow by a factor of e^^. Finally, if the moduli mass is of the 
same order as a typical mass taken in theories of quintessence, m ~ 10~^^Mp, we can have 
an accelerated expansion of the universe by a factor e^^, in a sub-Planckian regime of field 
space, just from trapping. (In this last case, as in ordinary quintessence models, tuning is 
required.) 

Thus, the stage of inflation in this simple model is shorter than the usual 60 e-folds, 
but it may nevertheless be very useful for initiating a first stage of inflation in theories 
where this would otherwise be impossible, or for diluting unwanted relics at the later stages 
of the evolution of the universe. Moreover, this scenario can easily describe the present 
stage of acceleration of the universe. 

One can also make the effect more substantial by constructing a more complicated 
scenario, consisting of a chain of particle production events at locations (p = (pi, where 
some fields Xi become light. The field (p may be trapped and enter the stage of parametric 
resonance near each of these points. Correspondingly, the universe enters the stage of 
inflation many times. One could arrange for 60 e-folds of inflation by taking, for example, 
m ~ 10^ GeV, ~ 10. 




(8.5) 
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Fig. 4: The D-brane picture of a series of trapping events. 

A D-brane example provides a useful geometrical model of this process. Suppose we 
have an observable brane B and another brane A approaching it. Suppose also that there 
are a number of other branes in between A and B. Each time the moving brane passes 
through one of the standing intermediate branes, stretched strings are created and slow 
the motion of A. The cumulative effect of a number of standing branes is perceived on the 
observable brane as a slowing-down of the motion of A due to the interactions. 

One should note that inflation in our scenario is rather unusual: the inflaton (p rolls a 
short distance, then oscillates for a long time, but with period much smaller than then 
rolls again, etc. This may lead to peculiar features in the spectrum of density perturbations. 
One can avoid these features if the points 4>i are very close to each other, and each of them 
does not stop the rolling of (f) but only slows it down. In this case, particle production will 
not lead to parametric resonance, so it is not very important to us whether the fields Xi 
are bosons or fermions, as long as their masses vanish at (pi. 

This scenario is similar to the string-inspired thermal inflation considered in [|12|] (see 
also [0), but our proposal does not require thermal equilibrium. The main effect which 
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supports inflation in our scenario is based on particle production and has a nonperturbative 
origin. (A closely-related mechanism uses the corrections to the kinetic terms in the strong 
coupling regime, where the particle production is suppressed 0.) We hope to return to a 
discussion of this possibility in a separate publication. 

9. Conclusion 

We have argued that the dynamics of rolling moduli is considerably modified due to 
quantum production of light fields. In fiat space quantum field theory, moduli typically 
become trapped in orbits around loci which have extra light degrees of freedom. In the 
presence of gravity, Hubble friction limits the field range the system samples, but any 
trapping events which do occur are enhanced by Hubble friction, which rapidly brings the 
modulus to rest at an ESP. Moduli trapping may aid in solving the cosmological moduli 
problem by driving moduli to sit at points of enhanced symmetry. Furthermore, the 
trapping of a scalar field which has a potential can lead to a short period of accelerated 
expansion in situations with steeper potentials than would otherwise allow this. Finally, 
the trapping effect has important consequences for the problem of vacuum selection, as 
it can reduce the problem to that of selecting one point within the class of ESPs. An 
intriguing feature of this process is that the trapping is more efficient near points with a 
large number of unbroken symmetries. 
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Appendix A. Particle Production Due to Motion on Moduli Space 

In this section we will calculate the quantum production of x particles, ignoring the 
effect of backreaction on the motion of (j). 

A mode of x with spatial momentum k obeys the wave equation 



[d^t +e + g''{ii'' + vh'')Yk = Q. (A.l) 



There are two solutions to this equation, u^j^ and w^'"*, associated to vacuum states with 
no particles in the far past and no particles in the far future, respectively. These two sets 
of modes are related by a Bogoliubov transformation 

ut = a^ul^' + f5kul^'\ (A.2) 

If we start in the state with no particles in the far past, then one can calculate the number 
density of particles in the far future to be 

rik = (A.3) 
in the /c*'* mode. This may be evaluated by solving equation (|A.1|) in terms of hypergeo- 



metric functions (see e.g. §3.5 of [|T3|), but we will present here a more physical argument. 

One can view ( |A.1|) as a one dimensional Schrodinger equation for particle scatter- 
ing/penetration through an inverted parabolic potential. If we send in a wave ij}]^ from 
the far right of the potential, part of it will penetrate to the far left, with an asymptotic 
amplitude Tjfi/'^"*, and part of it will be reflected back to the right, with an asymptotic 
amplitude Rkil^]^*, where and Rk are the transmission and reflection amplitudes.! 

The Bogoliubov coefficient in ( |A.2|) is determined in terms of these transmission and 
reflection amplitudes via 

Pk = ^. (A.4) 

Now we use a trick from quantum mechanics to relate R and T using the WKB 
method. If we are moving along the real time coordinate, the WKB form of the solution 
u)^{t) will be violated at small t, due to non-adiabaticity. However, if we take t to be 



The modes in the two problems are related by u]J^{t —>■ — oo) = T^il^l^** , etc. 
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complex then we can move from t = — oo to t = +00 along a complex contour in such a 
way that the WKB approximation 

~ ^ _ e-r v/^^T^^O^^T^^*' (A.5) 

is valid. Here the integral J* dt' becomes a contour integral along a semicircle of large 
radius in the lower complex t plane. For large \t\, we can estimate the phase integral in 
( |A.5| ) by expanding 

7 2 2 2 

Vfc^ + (72(^2 ^ ^2^2) r^gyt+ ^ . (A.6) 

^ y (/ C 

As we go around half of the circle, this term generates a factor 



^g-i7^^)-^(A;^+gV )/25t^-l/2 ^ ■^-7r{k^+g^fi^)/2gv_ f^^ ^) 

This is exactly the ratio between R* and T*, so we find 

nfc = = e--('='+9'M')/9-. (A.8) 

It is important to note that this result applies much more generally than for 

(f) = i^ + vt. (A. 9) 

In many cases the nonadiabaticity is only appreciable near the origin (p = 0, so that the 
near-origin trajectory can be approximated by (|A.9|) with some appropriate near-origin 
velocity v, even if the evolution away from the origin is very different from (p.2|) . 

Moreover, in analogous circumstances (T-dual in the brane context) with a nontrivial 
electric field, we obtain a similar expression due to Schwinger pair production; a related 
point was made in [^^. In addition, formula ( |A.8p applies not only to scalar fields, but 



also to fields of arbitrary spin. From this universal behavior, it is tempting to speculate 
that ( [A.8| ) could provide an effective model for string theory effects, but we will not pursue 
this direction here. 

The result ( [A.8| ) is nonperturbative in g (with g, not g'^, appearing in the denominator 
of the exponent); it is interesting to ask whether there is a simple interpretation of this 
nonanalytic, nonperturbative effect. Similarly, it is interesting to note that as discussed in 
§3, the potential for cp induced by particle production is linear, so that if extended to the 
origin it would have a nonanalytic cusp there. 
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Our results correspond to the low- velocity limit of the D-brane calculation by Bachas 
2^ . There, he obtains an imaginary part to the action for moving D3-branes of the form 



Im ^ oc V ^ ^ ( ^ ) exp(-n7^gfi^/v) (A.IO) 

^-^ n \Tvn/ 

n=l 

where we have translated his results into our variables. The first term in this expansion 
is proportional to the overlap / d^k\(3k\'^ giving the number density ( p.l4|) of produced 
particles; this agrees with what we expect from unitarity. More generally, backing away 



from this low-velocity limit, the calculation in combined with unitarity may provide a 
simple generalization of our results to the string case, as we briefly discussed in §5.2. 



Appendix B. Annihilation of the x Particles 

In this section we study the effects of collisions and direct decays of the created x 
particles, and demonstrate that for suitably chosen parameters the trapping effect receives 
only small corrections. More specifically, we place limits on the reduction of the x number 
density through processes like XX ~^ 4>4' ^md x ~^ VV^ where r] is some light field. 

Direct decays, if present, could easily ruin the trapping mechanism: if the x particles 
decay too rapidly into light fields then the energy stored in created x particles will not 
suffice to stop (/}. In this case the modulus will roll past the ESP, feel a transient tug 
toward the ESP while the x particles remain, and then gradually break free and glide off 
to infinity at a reduced speed. 

We will therefore consider only models in which couplings of the form xi^i^^ with 
•0, '0 very light, are negligible. As an example, one can easily exclude such decays in a 
supersymmetric model with a superpotential of the form W ~ g^X.^. Here X is a chiral 
superfield with scalar component x and fermion i/'^, and $ is a chiral superfield with scalar 
component (p and fermion t/jip. This generates Yukawa couplings of the form X'^x'^4> 
(pip^il}^, which do not allow decays from a component of X to purely $ particles. Thus, if 
all components of X are heavy, the X-particle energy density we produce cannot decrease 
by direct decays. In some of the simplest brane setups, exactly this situation is realized: a 
string which is heavy because it stretches between two branes separated in a purely closed 
string bulk space cannot decay perturbatively into two light, unstretched strings. 

2 

On the other hand, a priori we cannot ignore the coupling ^x?4>^ as it is this which 
gives rise to the desired trapping effect. This means that we must tolerate a certain rate 
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of annihilation (as opposed to direct decay). We will now review the cross section for this 
process and determine its effect on the number density appearing in ( p.lO| ). 



The Lorentz-invariant cross section for the annihilation process XX ~^ 4'4' is, written 
in terms of center of mass variables, 

where k and k' are the momenta of the ingoing and outgoing particles and E is the energy 
of the ingoing x particles. The reverse process (f)(f) — > xX tends to enhance the trapping 
effect. As we are in search of a lower bound on the number of x particles, we will simply 
omit this reverse process. 

We now determine the annihilation rate to find the rate at which x particles are lost. 
If we assume that all the x's are produced at t = 0, we find 



■h(k^ f\ r ^ ^ \ {kik2)'^-ml ^ ^ 

= - / dkMk2,t)^^—— -a{k^,k2). (B.2) 

n{ki,t) J -C/i-C/2 



Here w = y (A;i/c2)^ — m'^/EiE2 is the Lorentz-invariant relative velocity of the initial x's 
and cr(/ci, /C2) is the cross section, to be calculated using ( [B.l|) in the center of mass frame. 

We can simplify ( [B.2[ ) to get an upper bound on how fast x decays. Ignoring the 
momentum dependence on the right hand side of ( [B.2|) , which amounts to taking the 
non-relativistic limit, and ignoring the mass of 4> produced by the x particles, we have 

We can bound the integral in the second term on the right hand side by n^, the total 
number of x's produced, as given in ( p.l4|) . To approximate the time-dependence of the 



mass m^, we take = /x^ +v'^t'^, which is what the uncorrected motion for (j) would give. 
So we have finally 



which yields 



n{ki,t) 2n{^^ + vH^) 



nikA) ( g^n^ vt\ ,^ ^, 

> exp -- — ^ arctan — . (B.5) 
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This is clearly bounded from below by 



exp(-^). (B.6) 



SO that the number density is reduced over time by at worst the factor (|BT 
The total energy density in the x particles at a given time is therefore 



E= j dkn{k,t)y/k^ + g^iiJ^ + vH^) 

>n a/ g'^iiJ? + vH'^) exp f — — — — arctan — ) 
^ \ 2nnv /i / 



(B.7) 



From this we see that the mass amplification effect of the x particles inevitably prevails 
and stops (p from rolling arbitrarily far past the enhanced symmetry point. 

This reduction of the number density softens, but does not ruin, the trapping effect. 
Using the energy density ( [B.?] ) in the simple estimate leading to ( ^.12| ), we find a new 
estimate for 6^: 



4tc 



3 



g5/2 



Thus, although collisions never lead to an escape, they do lead to a somewhat increased 
stopping length cj)^. 

For suitably chosen parameters we can arrange that the effect of collisions is unimpor- 
tant and the estimates ( |2.12|) ,( p!8| ) approximately agree. For example, the final exponential 
factor, which encodes the consequences of annihilations, will be less important than the 
factor e^^^^ as long as g^v <^ /U^. 

We conclude that direct decays can be forbidden using symmetry, whereas collisions 
increase the stopping length ([B.8|) but do not ruin the trapping effect. 



Appendix C. Bending Angle and Energy Loss during the First Pass 

In this section we will demonstrate that one can obtain analytical results for the 
early stage of the motion of the moduli through a systematic expansion and iteration 
procedure. Specifically, we start with the uncorrected trajectory for the moduli, work out 
the production of x particles to the second order in w/u;^, feed this back into the equation 
of motion for the modulus, and then obtain the corrected trajectory for the modulus to 
leading nontrivial order of the iteration and Cj juP' . 
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The solution shows deflection and energy loss, as we would expect. However, without 
going to higher orders in the iteration, and without being able to detect particle production 
which is nonperturbative in cu/o;^ (as in (|2.9|) ), we do not expect our solution to capture the 
full dynamics. In particular, the leading-order calculations presented here do not capture 
the late-time behavior of the system, where any small particle production, such as that in 
(p^), can alter the trajectory from an escaping one to a trapped one. 

We choose to work with the simple model ( |2.1[) in which the complex modulus (f) = 
(pi + i(p2 has an ESP at the origin, where a real scalar fleld x becomes light. Without 
including backreaction, the trajectory ( p.2| ) leads to a time-varying frequency for the x 
excitations, 

ujkit) = ^/k-^ + g^iiJ^ + vH^). (C.l) 
We start by defining the generalized Bogoliubov coefficients and 



iiAt) = -i 



^a;(^) = TT^exp +2z / u;^^{t')dt' (3^^{t) 



(C.3) 



^kit) 

— 2 — 1 "fcV'J CT^P 1 / ^]^\'' I - fjj^^Lj cj'.y y-rt I ujj:yL JUL I 1 

(C.2) 

This is possible as long as uJk{t) is never zero. The consistency of this definition combined 
with the Klein-Gordon equation for leads to two coupled first-order ODE's for 
and f3j^{t) 

2u;^(t) 

The normalization for Uj^{t) implies 

alit) - Plit) = 1. (C.4) 
Finally, the initial condition for x{t) a,s t —>■ — oo corresponds to 

= 1, = 0. (C.5) 



Similar calculations appeared in M. 
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{x\t)) 



47r2 Jo ujk{t) 



When the frequency Ukif) changes adiabaticaUy, |ci;/a;^| ^ 1, the exponentials approxi- 
mate the positive and negative frequency modes of the system and [3k it) describes particle 
production. 

In terms of and I3^(t) 

l + 2|/3fc(t)|2 + 2Re(^afc(0/?:(0exp(^-2z j Ukit')^^^^ . 

(C.6) 

The first term in the above equation is the ordinary Coleman- Weinberg potential (after 
cancelling against some polynomial counter terms), and we ignore it for the reasons men- 
tioned in §3.2. The second and third terms encode the particle production effects; we will 
calculate these by an adiabatic approximation. 

To the leading nontrivial order, akif) and /3fc(t) satisfy 

=0 



where we have used the zeroth-order Bogoliubov coefficients, ct^it) = 1, = 0. These 

are subject to the initial conditions (|C.5|) , so that ak{t) = 1 and 



M) = ( ^ - ^ + ^ + - ) exp(-2z / u,dt'). {C.i 



This yields 



4ul AujI 8ul 



X {t)) - 6m = zr^ / dkk — — — — - - — — + 



2 2 

V n 



(C.9) 



487r2(^2 + t;2t2)2- 

To the first order of iteration, the corrected equation of motion for (p is 

2 2 2 

487r2(^2 _^ ^2^2)2^ V ^ 

with initial conditions 

(P{to) = vto + i^, ^{to)=v. (C.ll) 

We will eventually take the initial time to to be — oo. Moreover, we have in mind vq < 0, 
/U > 0, so that (p comes from (pi = +oo with a displacement in the positive direction of (p2 
and shoots towards the left. 
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The solution for the above system is 
where 



487r2 

^iW = ^^{z-t)'-^{z + t)'-^, (C.13) 

'^2(2) = 2 + 2 , 

and the branch cut goes from z = —i to z = i. 

As time goes from —00 to +00, z goes from 0"*" to 00 and then to 0~. We can enclose 
the contour by a semicircle either in the upper or lower half-plane, giving the same result: 



(C.14) 



We find that at late time 

(p{t —>■ +00) —>■ — vt(cos{n\^l + c) H 1= sin(7rVTT~c) ) 

V vl + c ^ 

— z/x^cos(7rvTT~c) + zVT+~csin(7rVTT~c) j (C.15) 

The leading large time behavior is captured by the first term, which describes a con- 
stant linear motion on the moduli space. It shows a deflection in the direction of the 



ESP (counterclockwise) determined by the coefficient of the leading term ( cos(7rvTT~c) + 



sin(7rvT+c) j . As we expect, if (7 = the modulus goes through a straight line 
from right to left. The bending angle increases monotonically, with no upper bound, as g 
increases. Thus, if g is very big, (p can orbit around the origin many times, with a growing 
radius, until eventually it shoots out to infinity in a fixed direction. 

The modulus loses a finite fraction of its energy as it passes by the ESP. This can 
be seen from the fact that the coefficient of —vt has norm less than one in (|C.15|) . It so 
happens that at this order of iteration the fraction of energy lost depends only on the 
coupling g: 

-AE/Eo = — ^ sin^rTrv^TTc). (C.16) 
1 + c 
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The above picture is valid only in certain parameter regimes. Specifically, if taking 
only the leading nontrivial contribution to (x^(t)) is to make any sense, then the non- 
adiabaticity parameter must be small: 



91^ 



< 1. 



(C.17) 



This condition also implies that the effects of the produced x particles on the motion of 
(/) shows up only at exponentially late times, so that it is sensible to have an intermediate 
stage in which the cj) field has a linear motion. Of course, we know this is not the full 
picture. The linear growth of the mass of the x particles would eventually stop cj) and 
bring it back. 



Appendix D. Classical Trapping Versus Quantum Trapping 

In this section we will compare our quantum trapping mechanism with the purely 
classical trapping proposed in 



Consider for simplicity a theory of two real scalar fields, cj) and x, with the interaction 

2 

^(p'^X^- III our discussion in the main text we assumed the initial conditions (x) = 0, {(p) ^ 
0, X = 0, = f . Potentially interesting classical dynamics arises in the more general case 
in which the initial velocity of x is nonzero |Q . 



Let us therefore consider the classical behavior of these fields, ignoring particle pro- 



duction entirely. If we define v = yx^ + 4''^, then energy conservation implies that the 
trajectory of (f) and x is bounded by the surface 9'^(px^ = . The fields will evidently 
start bouncing off the curved walls of the potential. This bouncing will be highly random. 
Naively, one would expect that on average the fields become confined in the region 

m = ix') - (D.l) 

This result would coincide with our estimate for the amplitude of the oscillations of (p at 
the end of the stage of parametric resonance, cf. (|3.21|) . 

However, the situation is more complicated. As we are going to show, the fields spend 
most of the time not at \(f)\ ~ |x| ~ but exponentially far away from this region, 

moving along one of the flat directions of the potential. 

To see this, note that because of the chaotic nature of the bouncing, it will occasionally 
happen that the flelds enter the valley x ^ ^ at a small angle to the flat direction, i.e. 
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with velocities obeying \x\ ^ 4>- Defining |x| = av, we are interested in the case that the 
angle a happens to be small. 

Energy conservation implies that the amplitude of the oscillations of x the initial 

2 

stage of this process is approximately Because of the interaction term ^^^x^; these 

2 2 

oscillations act on the field (f) with an average returning force ~ " " , which corresponds 



to the logarithmic potential V{(f)) ~ a^v"^ log^ Clearly, this potential will eventually 



pull the field (f) back to the ESP ^ = 0. However, this happens at exponentially large 
the field starts moving back only after its value approaches 



c/a 



e"/", (D.2) 

9 

where c = 0{1). 

Once again, because the bouncing process is highly random, we do not expect that 
the probability to enter the valley at a small angle a is exponentially suppressed. This 
means that after bouncing back and forth near the point = x = 0, the fields cf) and x 
eventually enter one of the valleys at a small angle, and subsequently spend a very long 
time there. In general, the fields will spend an exponentially long time at an exponentially 
large distance from the origin. Thus, the classical trapping mechanism, unlike the particle 
production mechanism described in our paper, does not lead to a permanent trapping of 
the fields in the vicinity of the point = x = 0. 
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